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ABSTRACT 

A  simple  formula  is  developed  in  this  paper  for  the  rela¬ 
tive  annual  damaging  ultraviolet  radiation  (DUV)  dose  at  dif¬ 
ferent  tropical  and  mid-latitude  sites.  The  formula  consists 
of  six  multiplicative  factors  which  include  the  effects  of 
amount  of  ozone,  latitude,  altitude,  cloudiness,  ground  albedo, 
and  amount  of  aerosols.  A  seasonal  ozone  variation  factor  is 
introduced  to  modify  the  tropical  relative  DUV  formula  for  ap¬ 
plication  to  mid-latitude  sites.  The  approach  involves  corre¬ 
lations  of  sometimes  sparse  data,  and  remains  to  be  validated 
in  the  general  sense. 

The  formula  should  be  useful,  where  more  exact  data  are 
not  available,  in  studies  of  the  effects  of  solar  ultraviolet 
radiation — and  its  possible  increase  from  a  reduction  in  strat¬ 
ospheric  ozone — on  land  and  marine  ecological  systems  and  skin 
cancer  incidence  in  white  Caucasian  populations. 
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SUMMARY 


A  formula  is  developed  in  this  paper  which  can  be  used  to 
compare  the  annual  damaging  ultraviolet  radiation  (DUV)  doses 
at  different  geographical  sites.  Relative  DUV  dose  plays  an 
important  role  in  studies  of  the  effects  of  solar  ultraviolet 
radiation — and  its  possible  enhancement  by  stratospheric  ozone 
depletion — on  marine  ecosystems,  agricultural  crops,  skin  can¬ 
cer  incidence  in  white  Caucasians,  etc.  Such  a  formula  also 
can  be  used  to  select  sites  for  future  solar  UV-B  measurements. 
The  formula  consists  of  six  multiplicative  factors  which  in¬ 
clude  the  effects  of  amount  of  ozone,  latitude,  altitude, 
cloudiness,  ground  albedo,  and  amount  of  aerosols. 

The  formula  is  based  on  the  assumption  that  D,  the  relative 
annual  damaging  ultraviolet  radiation  dose  (DUV)  at  a  given  site 
can  be  expressed  as  the  product  of  six  separable  multiplicative 
factors ,  i. e . , 


D  -  dt  Dl  Dh  Dc  Da  Dg  , 

where  the  subscripts  x,  L,  h,  C,  A,  and  $  refer,  respectively, 
to  the  average  amount  of  ozone,  latitude,  altitude,  average 
cloud  amount,  ground  albedo,  and  amount  of  aerosols.  The  refer¬ 
ence  value  of  D  is  unity,  corresponding  to  an  equatorial  sea- 
level  site  with  an  average  annual  amount  of  ozone  of  240  m. 
atm-cm,  no  clouds,  zero  ground  albedo,  and  standard  atmosphere. 
The  two  most  significant  factors,  D,.  and  DT ,  are  empirically 
shown  in  the  paper  to  be  independent  of  each  other.  The  remain¬ 
ing  factors  are  correction  factors  which  could  involve  some  in¬ 
terdependency  effects.  While  the  magnitudes  of  these  effects 
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are  probably  of  second  order,  the  validity  of  the  formula  re¬ 
mains  to  be  demonstrated  in  a  general  sense. 

The  expressions  derived  for  the  six  factors  are  as  follows: 
D  =  1  -  0.00484  (t  -  240) 

T 


in  the  tropics  and 

Dt  =  9.80  x  10-6  t2  -  1.0186  x  10-2  t2  +  2.886 

at  mid-latitudes  where  x  is  the  annual  average  ozone  thickness 
in  m.  atm-cm. 

D,  -  e’3'74  x  10'1' 

Li 

in  the  tropics  and 

D  .  e'3-74  x  10'11  ^  a  (t,  L) 

L 

at  mid-latitudes  where  L  is  the  latitude  in  degrees  and  a  (t,  L) 
is  a  correction  factor  (Eq.  30). 

D  =  1  +  0.060  h 
h 

where  h  is  the  altitude  in  degrees. 

Dc  =  1  -  0.50  C 

where  C  is  the  average  cloud  amount 

DA  =  1  +  0.50  A 
where  A  is  the  ground  albedo 

Dg  =  1  -  0.093  (8  -  1) 
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where  8  is  the  ratio  of  the  amount  of  aerosols  to  the  standard 
amount  of  aerosols. 

For  mi ’-latitudes ,  it  was  found  necessary  to  introduce  an 
additional  seasonal  ozone  correction  factor, 

1  +  p  (t,  L)  , 

where  p  (x,  L)  is  given  by  Eq.  (37)  on  p.  43  for  the  Southern 
Hemisphere  and  by  Eq.  (43)  on  p.  46  for  the  Northern  Hemisphere 

The  ozone  data  base  used  in  this  paper  was  the  Atlas  of  th 
Global  Distribution  of  Total  Ozone  July  1957 -June  1967  (Ref.  1) 
Calculations  of  DUV  doses  by  A.E.S.  Green  and  T.  Mo  in  the  CIAP 
program  (Refs.  2  and  11)  provided  basic  input  information  essen 
tial  to  many  of  the  empirical  relationships  derived. 

A  comparison  of  the  relative  annual  DUV  dose  was  made  for 
19  selected  tropical  sites,  considering  only  the  effects  of 
amount  of  ozone,  latitude,  and  altitude.  It  was  found  that 
Cerro  de  Pasco,  Peru,  a  mining  town  at  an  altitude  of  4.40  km 
ar.d  a  latitude  of  10° S,  and  Quito,  Ecuador,  on  the  equator  at 
an  altitude  of  2.85  km,  both  have  a  relative  annual  DUV  approx¬ 
imately  80  percent  higher  than  Key  West,  Florida,  which  is  a 
sea-level  site  at  a  latitude  of  24Js°N. 


I.  INTRODUCTION 


The  objective  of  this  paper  is  the  derivation  of  a  formula 
which  can  be  used  to  compare  annual  damaging  ultraviolet  radia¬ 
tion  (DUV)  dose  at  different  geographical  sites.  In  Chapter  II, 
a  formula  for  tropical  sites  is  derived,  and  a  modification  of 
this  formula  for  mid-latitude  sites  is  derived  in  Chapter  III. 

The  principal  factors  considered  are  amount  of  ozone  (x),  lati¬ 
tude  (L),  altitude  (h),  cloudiness  (C),  ground  albedo  (A),  and 
amount  of  aerosols  (8).  In  this  exploratory  effort,  it  was 
found  that  for  tropical  sites,  these  six  factors,  while  not  all 
independent  of  one  another,  can  to  a  first  approximation  be  com¬ 
bined  in  a  product  formula,  i.e., 

D  -  dt  dl  Dh  Dc  da  ce  •  (1) 

where  D  is  the  relative  annual  damaging  ultraviolet  radiation 
dose  and  each  of  the  component  factors  are  referenced  to  an  ap¬ 
propriate  baseline  value  of  unity.  All  relative  DUV  doses  used 
in  this  paper  represent  the  solar  energy  incident  on  a  horizontal 
surface,  weighted  by  the  erythemal  response  spectrum. 

A  formula  for  mid-latitude  sites  was  derived  by  modifying 
two  of  the  above  six  multiplicative  factors  and  introducing  an 
additional  seasonal  ozone  variation  factor  (Chapter  III). 

Tropical  sites  are  of  special  interest  in  considering  the 
potentially  harmful  effects  of  a  reduction  in  stratospheric 
ozone  induced  by  fluorocarbons,  future  stratospheric  aircraft 
fleets,  etc.  Near  the  equator  the  protective  ozone  thickness 
has  a  relatively  stable  minimum  value  and  the  sun  a  high 
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noontime  elevation  angle  throughout  the  year.  Consequently,  a 
large  portion  of  the  earth's  surface  receives  a  maximal  annual 
dose  of  biologically  damaging  ultraviolet  (DUV)  radiation  (42 
percent  of  the  earth's  surface  lies  between  latitudes  25°N  and 
25° S ) .  In  the  event  of  a  significant  reduction  in  global  strat¬ 
ospheric  ozone,  the  equatorial  zone  could  be  confronted  with 
levels  of  ultraviolet  radiation  perhaps  not  experienced  on  earth 
for  billions  of  years.  While  it  is  possible  to  equate  the  fu¬ 
ture  new  level  of  ultraviolet  radiation  at  a  sea-level  mid-lati¬ 
tude  region  to  that  currently  received  at  a  sea-level  latitude 
closer  to  the  equator,  this  latitude  translation  cannot  be  made 
at  the  equator.  Migration  for  stressed  equatorial  ecological 
systems  at  sea  level  is  therefore  not  a  possibility.*  It  would 
be  highly  desirable  to  investigate  the  behavior  of  ecological 
systems  in  those  equatorial  regions  of  maximal  ultraviolet  radi¬ 
ation.  The  initial  motivation  for  this  paper  was  the  derivation 
of  a  simple  formula  which  can  be  used  to  identify  those  low- 
latitude  sites  which  today  are  receiving  the  highest  annual 
damaging  ultraviolet  radiation  doses  on  earth. 

While  tropical  sites  are  of  great  ecological  interest,  the 
vast  majority  of  the  white  Caucasian  populations,  which  are  sus¬ 
ceptible  to  solar-ultraviolet-induced  skin  cancer,  are  to  be 
found  living  at  mid-latitude  sites.  A  quick  method  of  calculat¬ 
ing  annual  relative  DUV  dose  may  be  of  some  use  in  comparing 
skin  cancer  incidence  at  various  geographic  sites.  However,  it 
must  be  recognized  that  while  many  ecological  systems  can  be  ex¬ 
posed  daily  to  the  full  DUV  dose  on  the  ground,  humans  receive 
only  a  fraction  of  the  dose,  depending  on  individual  exposure 
and  clothing  habits. 

if - 

Downward  altitude  migration  of  ecosystems  that  have  adapted  to 
higher  elevations  could  provide  some  limited  relief.  Marine 
animals  might,  where  depth  permits,  move  to  deeper  waters. 


II.  TROPICAL  SITES 


A.  OZONE 

The  thickness  of  the  ozone  column,  x,  varies  with  latitude. 
In.  Pig.  1  (from  Ref.  1)  is  shown  the  global  distribution  of 
total  ozone  averaged  over  the  10-year  period  starting  in  July 
1957.  The  worldwide  distribution  of  the  stations  used  to  ob¬ 
tain  data  is  also  indicated  in  the  figure.  The  contour  map  in¬ 
dicates  that  the  northern  part  of  South  America  and  the  central 
part  of  Africa  had  the  minimum  average  total  ozone  column  of 
approximately  240  x  1CT3  cm,  or  240  m.  atm-cm,  for  the  period 
1957-1967.  If  only  the  ozone  factor  is  considered,  the  highest 
ultraviolet  radiation  levels  could  be  expected  to  be  found  with¬ 
in  the  240  m.  atm-cm  contour  of  Fig.  1. 

The  use  of  an  average  value  of  total  ozone  in  the  formula 
derived  below  for  the  tropics  would  be  inappropriate  for  regions 
outside  the  tropical  zone  because  of  the  large  ozone  fluctuations 
with  season.  The  effect  of  season  on  amount  of  ozone  is  minimal 
near  the  equator,  as  indicated  in  Fig.  2  (from  Ref.  1).  Also, 
differences  in  total  average  ozone  of  only  approximately  10  m. 
atm-cm  are  to  be  found  at  the  equator,  while  a  difference  of 
70  m.  atm-cm  is  found  between  longitudes  130°E  and  40°E  at  a 
latitude  of  50°N  (Fig.  1). 

The  ozone  function  Dx  at  the  equator  is  shown  as  a  function 
of  x  in  Fig.  3.  The  five  circled  points  shown  are  based  on  the 
sum  of  the  12  monthly  DUV  tabulated  values*  calculated  by  T.  Mo 
and  A.E.S.  Green  (Ref.  2)  at  sea  level  for  a  standard  amount  of 


Includes  both  scattered  and  direct  solar  radiation  incident 
on  a  flat  horizontal  surface. 
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FIGURE  1.  Global  distribution  of  total  ozone 
averaged  over  the  period  1957-1967. 
(Source:  Ref.  1) 


N  LATITUDE  S 


FIGURE  2.  Latitudinal  variation  of  aver¬ 
age  total  ozone  averaged  over 
the  period  1957-1967.  (Source 
Ref.  1) 
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FIGURE  3.  Relative  DUV  vs  ozone  thickness  at  the  equator 
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aerosols  in  a  clear  atmosphere  with  zero  ground  albedo.  A  para¬ 
bolic  fit  to  the  three  calculated  points  at  256,  288,  and  320  m. 
atm-cm  was  made,  leading  to 

Dx  -  9.^24  x  10~H  x2  -  0.830  T  +  206.2  .  (2) 

However,  a  simpler  linear  approximation  can  be  used  here  to 
cover  the  narrow  band  of  expected  values  of  t  near  the  equator. 
The  following  linear  equation  is  adequate  for  the  purpose  of 
this  paper  to  cover  the  tropical  range  of  240  -  280  m.  atm-cm 
(10-year  average). 


Dt  =  62  -  0.30  (t  -  240)  .  (3) 

Dividing  by  62,  we  have  the  relative  ozone  value  function 

D  -  1  -  0.00484  (t  -  240)  .  (4) 

Eq.  (4)  can  be  used  to  compare  the  annual  DUV  at  two  sites  of 
different  longitude  along  the  equator  (L  =  0°)  if  they  differ 
in  their  t  value  and  all  of  the  other  four  factors  are  assumed 
to  be  equal. 

The  second  term  in  Eq .  (4)  gives  the  fractional  decrease 
in  moving  from  a  site  having  240  m.  atm-cm  of  ozone  to  one  with 
a  higher  value  of  x.  For  every  1  percent  increase  in  ozone 
along  the  equator,  there  is  a  1.16*  percent  decrease  in  DUV  (as 
compared  to  approximately  2  percent  decrease  for  mid-latitude 
sites) . 

B.  LATITUDE 

In  Fig.  4  the  variation  of  relative  DUV  with  the  square  of 

p 

the  latitude  (L  )  is  plotted  on  semi-log  paper  for  values  of  x 

In  agreement  with  Fig.  4.4,81  in  The  Report  of  the  Committee 
on  Meteorological  Effects  of  Stratospheric  Aircraft  (COMESA) 
1972-1975,  Part  2,  1975. 
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between  256  and  320  m.  atm-cm.  The  circled  points  are  again 

based  on  the  calculations  of  T.E.  Mo  and  A.E.S.  Green  (Ref.  2). 

There  are  two  significant  empirical  observations  to  be  made  in 

Fig.  4:  (1)  for  a  given  value  of  t  the  relative  DUV  values  fall 

almost  exactly  on  a  straight  line  for  latitudes  less  than  25°, 

and  (2)  the  slopes  of  the  lines  are  almost  exactly  equal.  With 

these  two  fortuitous  observations,  it  is  possible  to  accurately 

formulate  the  relative  annual  DUV,  considering  only  ozone  and 

latitude,  as  the  product  D  DT  where 

T  Li 

Dl  .  e’3-711  x  10'“  L2  ,  (5) 

and  L  is  in  degrees. 

With  Eqs.  (4)  and  (5)  it  is  possible  to  compare  the  rela¬ 
tive  DUV  for  any  two  tropical  sites  in  Fig.  1,  assuming  the 
other  four  factors  are  equal.  Thus,  for  example,  the  coastal 
town  of  Townsville  in  Queensland,  Australia  at  a  latitude  of 
19°S  with  a  x  of  260  m.  atm-cm  (Fig.  1)  had,  over  the  period 
1957-1967,  an  average  DUV,  relative  to  an  equatorial  sea-level 
site  in  South  America  or  Africa,  of 

(1  -  0.0968)  e"0,135  =  0.78  >  (6) 

assuming  equality  of  the  other  three  factors,  i.e.,  cloudiness, 
ground  albedo,  and  amount  of  aerosols. 

C.  ALTITUDE 

The  influence  of  site  altitude  on  the  total  (direct  +  dif¬ 
fuse)  ultraviolet  flux  received  at  the  ground  has  been  investi¬ 
gated  by  S.V.  Venkateswaran  et  al.,  (Ref.  3).  In  Fig.  5  is 
shown  the  flux  normalized  to  sea  level  for  wavelengths  of  300.4 
and  311.4  nm  for  different  solar  zenith  angles  and  for  altitudes 
up  to  2  km  (Ref.  3).  As  indicated  in  Fig.  5,  flux  was  found  to 
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FLUX  NORMALIZED  TO  SEA  LEVEL 


FIGURE  5.  The  influence  of  the  altitude  of  the  station  on 
the  total  flux  received  at  the  ground  is  illustrated  for 
two  wavelengths  (300.4  and  311.4  nm)  and  for  different 

solar  zenith  angles  0Q  (=  cos"1  yQ).  The  fluxes  are  nor 

malized  to  sea-level  values.  The  ground  albedo  is  as¬ 
sumed  to  be  zero.  (Source:  Ref.  3) 


FIGURE  6.  Ratio  of  flux  at  1  km  altitude  to  flux  at  sea 
level  vs  solar  zenith  angle 


increase  almost  linearly  with  altitude,  depending  on  wavelength 
and  solar  zenith  angle.  While  ground  albedo  was  assumed  to  be 
zero  in  these  calculations,  it  was  stated  that  the  dependence 
of  the  ground  albedo  effect  on  site  altitude  is  small  for  mod¬ 
erate  site  elevations. 

The  relative  DUV  dose  D  (h,  6)  is  a  function  of  altitude  h 
and  solar  zenith  angle  0.  To  obtain  a  value  independent  of  0 
it  is  necessary  to  integrate  over  0.  At  the  equator  the  solar 
zenith  angle  0  is  close  to  being  uniformly  distributed*  between 
0°  and  90°  because  of  the  constant  rotation  of  the  earth. 
Therefore,  the  relative  altitude  DUV  function  can  be  approxi¬ 
mately  calculated  from  the  equation 


'h 


*  1  +  h 


n/2 

/ 


/ 


D ( 0 ,  0) 

T72 


[d(1,  0)  ,1 

[dTo^TT 


d0 


(7) 


D(0,  0)  d0 


where  D(l,  0)  is  the  relative  DUV  at  an  altitude  of  1  km  and 
for  solar  zenith  angle  0.  The  function  is  plotted  vs  0 

in  Fig.  6  using  the  information  in  Fig.  5  for  X  =  300.4  nm  and 
X  =  311.4  nm.  Neither  of  these  wavelengths  represents  the  DUV, 
but  a  wavelength  of  approximately  305  nm  lying  midway  between 
the  curves  would  be  representative  (Ref.  4). 

In  Fig.  7  the  relative  DUV  function  D(0,  0)  is  plotted  as 
a  function  of  the  solar  zenith  angle  for  amounts  of  ozone  be¬ 
tween  240  and  400  m.  atm-cm  (Ref.  5).  When  plotted  on  semi-log 
paper  vs  O2*'*  (Fig.  8),  it  is  found  that  D(0,  0)  is  excellently 
approximated  by 


0(0,  0) 


e-7 .2355  x  10"5  02,5 


(8) 


for  0  <  60°  and  240  <  t  <  400  m.  atm-cm. 

1 - 

Exactly  uniformly  distributed  on  two  days  of  the  year  at  the 
equator. 


ctiveness  on 
or  five  ozone 


A  check  on  the  validity  of  Eq.  (8)  is  shown  in  Pig.  9.  The 
curve  of  relative  DUV  vs  time  of  day,  as  measured  by  a  Robert¬ 
son  Berger  meter  on  a  clear  day  at  Mauna  Loa,  Hawaii  on  August 
11,  1974  (Ref.  6),  is  seen  to  match  remarkably  closely  the  Eq. 
(8)  values  or  the  empirical  curve  of  Fig.  7.  To  locate  the  Eq. 
(8)  values  in  Fig.  9,  the  maximum  solar  zenith  angle  was  esti¬ 
mated  to  have  occurred  at  12:12  PM  and  each  10°  of  solar  zenith 
angle  was  assumed  to  correspond  to  43.5  minutes  of  time.* 

The  Mauna  Loa  observatory  is  located  on  a  mountain  top  at 
an  elevation  of  approximately  4  km.  Fig.  9  suggests  that  the 
dependence  of  relative  DUV  on  solar  zenith  angle  8  is,  to  a 
first  approximation,  independent  of  altitude  h,  i.e., 

D(h ,  6)  *  k  D(0,  8)  (9) 

where  k  is  a  constant.  Numerical  integration  of  Eq.  (7),  using 
Fig.  6  and  the  X  -  300.4  nm  curve  of  Fig.  5,  yielded 

Dh  *  1  +  0.0605  h  .  (10) 

Since  the  X  =  311.4  nm  curve  is  approximately  1  percent  below 
the  X  =  300.4  nm  curve  in  Fig.  6,  the  interpolated  Dh  function 
for  the  DUV  curve  at  the  equator  is  estimated  to  be 

Dh  =  1  +  0.055  h  .  (11) 

In  view  of  the  approximation  in  Eq.  (7)  that  assumed  the 
sun  reaches  the  zenith  every  day,  and  that  furthermore  an  equa¬ 
tion  is  desired  which  is  representative  of  non-equatorial  lati¬ 
tudes,  Eq.  (11)  is  rounded  off  in  an  upward  direction  to 

Dh  *  1  +  0.060  h  .  (12) 

1 - 

The  length  of  day  at  20°N  latitude  is  13.03  hr  on  August  11. 
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An  error  in  the  coefficient  of  h  would  have  a  maximum  effect  on 
the  value  of  at  a  maximum  value  of  h.  Hence,  at  h  =  4  km,  a 

5  percent  error  in  the  coefficient  of  h  would  correspond  to  only 
a  1  percent  error  in  the  value  of  D^.  It  follows  that  a  good 
approximation  for  the  relative  DUV  dose  for  sites  at  latitudes 
below  25°  is  given  by 


D„  Dt  D.  =  [1  -  0.00484  (t-240)] 

T  "  “  7ii  x 

(1  +  0.060  h)  e  5,1 


(13) 


assuming  cloudiness,  ground  albedo,  and  amount  of  aerosols  are 
equal . 


Equation  (13)  may  be  used  to  find  the  combination  t,  L,  and 
h  values  that  will  match  the  unity  value  of  relative  DUV  at  an 
equatorial  site  (L  =  0°)  at  sea  level  (h  =  0)  with  amount  of 
ozone  t  =  240  m.  atm-cm.  The  site  altitude  h  required  is  illus¬ 
trated  in  Fig.  10  as  a  function  of  latitude  L  with  t  as  parameter. 
Also  sketched  in  Fig.  10  is  the  approximate  altitude  limit  for 
biological  ecosystems  of  approximately  5  km,  and  the  approxi¬ 
mate  upper  and  lower  limits  of  ozone  amount  vs  latitude  as  given 
by  Fig.  1.  If  L  =  15°  and  t  =  255  m.  atm-cm,  the  altitude  re¬ 
quired  to  match  the  relative  DUV  of  a  sea-level  equatorial  site 
with  t  =  240  m.  atm-cm,  is  approximately  3  km.  Altitudes  con¬ 
siderably  higher  than  3  km  are  readily  found  in  populated  re¬ 
gions  of  Peru  and  Bolivia. 


D.  CLOUDINESS 

The  effect  of  site  cloudiness  on  relative  DUV  is  very  com¬ 
plex  but  it  is  possible,  as  shown  below,  to  use  some  available 
data  to  derive  a  simple  approximate  formula  for  D~. 

In  Table  1  is  listed  the  annual  DUV  count,  as  measured  by 
the  Robertson-Berger  meter  for  10  U.S.  sites  for  the  year  1974 
(Ref.  6).  The  Minneapolis  daily  DUV  count  for  1974  is  illus¬ 
trated  in  Fig.  11  (Ref.  6).  The  complex  effect  of  clouds  is 
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FIGURE  10. 


Minnesota 


demonstrated  by  the  wild  deviations  from  the  envelope  of  the 
sine-like  curve.*  On  August  3  the  DUV  dose  was  one-eighth  of 
the  clear-day  dose  (Ref.  6).  The  days  when  the  DUV  count  fell 
on  the  envelope  of  the  curve  were  clear  days.  By  integrating 
the  area  of  the  curve  under  the  envelope  it  is  possible  to  cal¬ 
culate  the  value  the  Robertson-Berger  meter  would  have  measured 
had  there  been  no  clouds  over  the  site  for  all  365  days  in  197 ^ . 
This  integral  is  also  tabulated  in  Table  1  as  calculated  from 
the  curves,  similar  to  Fig.  11,  which  are  available  for  each  of 
the  10  sites  in  Ref.  6.  The  ratio  of  the  measured  annual  DUV 
value  to  the  integral  yields  the  desired  cloudiness  factor  Dc. 
Note  that  in  197 4  the  value  of  Dc  ranged  from  a  low  of  0.67  in 
Philadelphia  to  a  high  of  0.90  in  Albuquerque. 

TABLE  1.  CLOUD  FACTOR  FOR  10  U.S.  SITES  IN  1974 

ANNUAL  DUV  f  * 

AREA  LATITUDE  (°N)  COUNT  x  10~6  J  DC 


Mauna  Loa,  HI 

19.5 

2.77 

3.24 

0.85 

Tallahassee,  FL 

30.4 

1.66 

2.16 

0.77 

El  Paso,  TX 

31.8 

2.24 

2.52 

0.89 

Fort  Worth,  TX 

32.8 

1.61 

2.20 

0.73 

Albuquerque,  NM 

35.1 

1.89 

2.10 

0.90 

Oakland,  CA 

37.7 

1.51 

1  .83 

0.83 

Philadelphia,  PA 

39.9 

1.11 

1.65 

0.67 

Des  Moines,  IA 

41.5 

1.25 

1.75 

0.72 

Minneapolis,  MN 

44.9 

1.07 

1.53 

0.70 

Btsrmarcky  NO**  *  *  ■ 

*45.8  -  *  * 

'*‘1.13*  '  *  ' 

1  1 64 

0.69 

* 

Annual  DUV  count 

with  no  clouds. 

ft 

Only  a  very  small 

fraction  of  the 

deviations  can 

be  attributed 

to  dally  ozone  thickness  variations. 

17 


A  comparison  of  frequency  of  cloudiness  by  cloud  amount  for 
the  year  1974,  the  same  year  as  that  represented  by  the  data  in 
Table  1,  is  available  for  two  of  the  sites  in  Table  1  (Fort  Worth 
and  Minneapolis)  from  calculations  made  by  R.  Penndorf  (Ref.  7) 
and  is  reproduced  in  Table  2. 


TABLE  2.  COMPARISON  OF  CLOUDINESS  IN 
FORT  WORTH  AND  MINNEAPOLIS 


CLOUD  AMOUNT 


PERCENTAGE  OF  OCCURRENCES 


0.0  -  0.1 

38 

30 

0.2  -  0.8 

16 

10 

0.9  -  1.0 

46 

60 

A  linear  relationship  of  the  relative  DUV  factor  with 
cloudiness  was  found  by  K.  Buttner  in  1938  (Ref.  8).  As  indi¬ 
cated  in  Fig.  12,  the  cloudiness  factor  Dc  can  be  well  approxi¬ 
mated  by  the  linear  equation 


Dc  -  1  -  0.55  C  ,  (14) 

where  C  is  the  average  cloud  amount  (unity  for  a  complete 
overcast ) . 

With  the  more  accurate  data  of  Tables  1  and  2  (there  were 
no  DUV  meters  in  1938)  it  is  now  possible  to  check  on  the  valid¬ 
ity  of  Eq.  (14). 

If  is  defined  as  the  average  cloud  amount  for  the  range 
of  cloud  amount  0.0  -  0.1,  i.e.,  0.05*,  C2  for  the  range  0.2  - 
0.8,  i.e.,  0.5,  and  for  the  range  0.9  -  1.0,  i.e. ,0.95,  then 
the  average  cloud  amount  for  the  year  may  be  represented  by 

1=3 

C  "  £  Ci  Pi  (15) 

1-1 

T - — 

Assumes  a  uniform  probability  distribution  over  the  cloud 
amount  interval. 


18 


FIGURE  12.  Dependence  of  relative  intensity 
of  UV-B  total  solar  radiation  on 
cloudiness.  (Source:  Ref.  8) 


where  pj^  is  the  frequency  of  occurrence  associated  with  average 
cloud  amount  as  given  by  the  values  in  Table  2. 

If  y  is  defined  to  be  the  coefficient  of  C  in  Eq.  (14), 
then  the  relative  cloudiness  factor  for  Fort  Worth,  using  the 
data  in  Tables  1  and  2  and  Eq.  (15),  is  given  by 

0.73  =  1  -  y  [0.38  x  0.05  +  0.16  x  0.5  +  0.46  x  0.95]  (16) 

and  for  Minneapolis  by 

0.70  =  1  -  Y  [0.30  x  0.05  +  0.10  x  0.5  +  0.60  x  0.95]  .  (17) 
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Solutions  of  Eqs.  16  and  17  yield  y  *  0.50  for  Port  Worth  and 
y  ■  0.47  for  Minneapolis.  Both  values  are  lower  than  the  0.55 
value  of  Buttner,  but  a  value  of  0.50  is  in  good  agreement  with 
the  straight  line  that  Buttner  could  have  drawn  to  better  fit 
the  points  shown  in  Pig.  9  (note  that  two  of  his  data  points  lie 
below  his  line  and  five  lie  above).  Thus,  the  cloudiness  factor 
adopted  here  is  given  simply  by 

Dq  =  1  -  0.50  C  .  (18) 

A  10  percent  error  of  0.05  in  the  coefficient  y  would  result  in 
an  error  of  only  1.5  percent  in  Dc  for  a  typical  C  value  of  0.25. 
Since  Eq.  (18)  was  derived  from  data  at  mid-latitude  sites,  it 
is  possible  that  a  slightly  different  coefficient  for  C  would  be 
found  for  tropical  sites. 

There  is  a  practical  difficulty  in  applying  Eq.  (18)  to 
compare  relative  DUV  of  low-latitude  sites.  Cloudiness  data  in 
the  form  given  in  Table  2  for  Minneapolis  and  Fort  Worth  are  not 
readily  available.  Photographs  taken  by  meteorological  satel¬ 
lites  have  been  used  to  compile  global  charts  of  relative  cloud 
cover  (Ref.  9).  These  charts  were  considered  not  suitable  for 
use  in  this  paper  because  of  the  inability  of  the  vidicon  camera 
to  distinguish  cloud  cover  from  snow  cover,  ice,  and  desert  ter¬ 
rain,  as  well  as  inadequate  resolution  (40  km).  Other  global 
charts  of  mean  cloudiness  (e.g..  Handbook  of  Meteorology  by 
Berry,  Bollay,  and  Beers,  McGraw-Hill,  1945)  are  too  crude  for 
the  purposes  of  this  paper. 

There  exist  other  intractable  meteorological  factors  that 
may  be  peculiar  to  a  given  site  and  significantly  influence  DUV. 
An  example  is  the  frequent  morning  fog  over  San  Francisco.  Ex¬ 
perimental  data  on  the  effect  of  fog  on  relative  DUV  are  not 
presently  available. 


E.  GROUND  ALBEDO 


Another  complex  factor  affecting  relative  DUV  is  the  ground 
albedo.  A  site  surrounded  by  snow-capped  mountains  will  receive 
more  ultraviolet  radiation  than  one  not  surrounded  by  snowy 
terrain  due  to  the  high  reflectivity  of  snow  (Ref.  4).  A  snow 
surface  may  have  a  reflectivity  ranging  from  0.4  to  0.6  in  the 
near  UV  (Ref.  3).  Reflection  of  ultraviolet  radiation  from 
vegetation  varies  with  the  type  of  plant.  The  UV  reflectance 
of  a  citrus  upper-leaf  surface  is  4  percent,  whereas  the  reflec¬ 
tance  of  the  upper-leaf  surface  of  a  century  plant  is  6  percent 
(Ref.  10).  Coastal  sites  will  have  different  ground  albedos 
than  interior  sites.  Albedos  for  water  surfaces  are  expected 
to  be  below  0.1  (Ref.  3).  City  sites  will  have  different  ground 
albedos  than  rural  sites.  Unfortunately,  measurements  of  ground 
albedo  in  the  ultraviolet  for  various  types  of  terrain  in  equa¬ 
torial  regions  (or  for  that  matter  almost  anywhere  else)  are 
unavailable.  Nonetheless,  the  effect  of  ground  albedo  on  ultra¬ 
violet  flux  can  be  and  has  been  investigated  (Ref.  3). 

In  Fig.  13  is  shown  the  effect  of  ground  albedo  on  total 
flux  for  three  ultraviolet  wavelengths  (300.4,  311.4,  and  339*8 
nm)  and  three  different  values  of  ground  albedo  (0.2,  0.4,  and 
0.6)  as  function  of  ozone  amount  for  a  sea-level  site. 

In  Fig.  14  the  effect  of  ground  albedo.  A,  on  relative  DUV 
is  shown  for  the  wavelengths  and  ozone  amounts  of  interest  in 

the  tropics.  Evidently,  the  narrowness  of  the  bands  of  interest 

results  in  a  dependency  which  is  essentially  independent  of  X 
and  t  can  be  well  represented  by  the  relationship 

DA  =  1  +  0.50  A  (19) 

.  '*  "  ■  -  *  ■»  •  •• 

for  expected  average  values  of  A  <  0.3.  It  is  a  curious  coinci¬ 

dence  that  in  Eq.  (19)  the  ground  albedo.  A,  has  the  same  coef¬ 
ficient  as  the  average  cloud  amount,  C,  in  Eq.  ( 1 8 ) ,  except  for 
sign. 


1.40 


X«S0Q4r>m 
X*3l  l.4nm 
X-339.8nm 


FIGURE  13.  The  effect  of  ground  albedo  (A)  on 

the  total  flux  (Fy)  for  three  differ¬ 
ent  wavelengths.  The  fluxes  for  three 
albedos  are  shown  normalized  to  the 
zero  albedo  values,  as  the  ozone  amount 
(Si)  over  a  sea-level  station  is  varied. 
(Source:  Ref.  3) 


FIGURE  14.  Effect  of  ground  albedo  on  ultraviolet 
flux 
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Because  of  the  flatness  of  the  curves  In  Pig.  13  for  A  = 

0.2,  it  can  be  seen  that  Eq.  (19)  is  also  a  fairly  good  approxi¬ 
mation  for  mid-latitude  sites  where  t  values  exceed  280  m.  atm-cm. 

F.  AEROSOLS 

The  effect  of  aerosols  on  relative  DUV  suffers  from  the 
same  complexity  difficulties  as  cloudiness  and  ground  albedo, 
depending  on  the  shape,  size,  and  size  distribution  of  the  par¬ 
ticles,  their  altitude  distribution,  and  their  index  of  refrac¬ 
tion,  as  well  as  the  obviously  important  particle  density 
(Ref.  3).  The  effect  of  aerosols  on  the  radiation  field  can 
therefore  only  be  calculated  for  particular  models. 

The  model  used  by  A.E.S.  Green  and  T.  Mo  (Ref.  11)  involved 
a  standard  amount  of  aerosols,  on  which  their  calculations  were 
based.  Calculations  were  made  for  0,  1,  2,  and  4  times  the 
standard  aerosols  as  indicated  in  the  DUV  curves  of  Fig.  15  for 
a  latitude  of  30°.  Note  that  the  DUV  curves  vary  essentially 
linearly  with  amount  of  aerosols  between  the  standard  amount 
(1)  curve  and  the  4-times-standard-amount  curve.  It  follows 
that  the  effect  of  aerosols  on  relative  DUV  can  therefore  be 
represented  approximately  by  the  simple  expression 

Dg  =  1  -  0.093  (B  -  1),  (20) 

where  8  is  the  ratio  of  the  amount  of  aerosols  at  a  site  to  the 
standard  amount  of  aerosols.  Equation  (20)  is  normalized  to  a 
value  of  unity  for  8=1.  A  10  percent  increase  in  the  amount 
of  aerosols  above  the  standard  amount  would  therefore  correspond 
to  a  1  percent  reduction  in  the  relative  DUV,  assuming  the  other 
five  factors  were  held  constant.  Since  relative  differences  of 
amount  of  aerosols  over  near-equatorial  sites,  where  industrial 
sources  of  atmospheric  pollution  are  relatively  rare,  can  be 
expected  to  be  small,  it  seems  highly  likely  that  the  aerosol 
relative  DUV  factor  is  the  least  significant  of  the  six  factors 
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FIGURE  15.  Dally  erythema  doses  with  four  assumed  amounts 
of  aerosols.  Amounts  of  aerosols  are  0,  1,  2, 
and  4  times  the  standard  aerosols  and  are 
labeled  on  the  curves.  One  MED  (minimum  7 
erythemal  dose)  is  assumed  to  be  200  (J/nr)  . 
(Source:  Ref.  11)  e 

considered  for  the  tropics.  This  is  supported  by  the  calcula¬ 
tions  for  four  U.S.  cities  (Fig.  16)  comparing  DUV  with  avail¬ 
able  turbidity  data  and  standard  aerosols  (Ref.  2).  Note  that 
the  effect  of  aerosols  at  mid-latitude  sites  is  greatest  in  the 
summer  months,  and  is  minor  for  the  small  cities  of  Boulder, 
Colorado,  and  Tucson,  Arizona,  compared  to  the  large  cities  of 
Los  Angeles  and  Chicago. 
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FIGURE  16.  Calculation  of  the  UV  monthly  erythema 
doses  for  four  metropolitan  areas  with 
the  data  of  local  turbidity  incorporated. 
Also  shown  are  the  results  obtained  with 
standard  aerosols  (solid  curves).  One 
MED  is  assumed  to  be  200  (J/m*) 

(Source:  Ref.  2)  e 


G.  FORMULA  FOR  RELATIVE  ANNUAL  DUV  AT  TROPICAL  SITES 

Substituting  Eqs .  13,  18,  19,  and  20  in  Eq.  (1)  gives,  for 
the  relative  annual  DUV  at  tropical  sites,  the  formula 

_i|  2 

D  =  [1  -  0.00484  (t  -  240)]  (1  +  0.06  h)  e'3*74  X  10  L 


x  (1  -  0.50  C)(l  +  0.50  A)  [1  -  0.093  (6  -  1)] 


(21) 


where 


x  =  amount  of  ozone  in  m.  atm-cm 
h  =  altitude  of  site  in  km 
L  =  latitude  of  site  in  degrees 

C  =  average  cloud  amount  (unity  for  complete  overcast) 

A  =  ground  albedo 

8  =  ratio  of  amount  of  aerosols  to  standard  amount  of 
aerosols . 

Eq.  (21)  can  be  used  without  much  loss  of  accuracy  for  sub¬ 
tropical  sites  but  for  latitudes  greater  than  30°  the  more  com¬ 
plex  mid-latitude  formula  Eq.  (44)  is  required. 

H.  RELATIVE  ANNUAL  DUV  AT  TROPICAL  SITES 

A  comparison  of  calculated  relative  annual  DUV  is  given  in 
Table  3  for  19  selected  tropical  sites.  Only  the  factors  for 
which  appropriate  data  is  available  were  included  in  this  compari¬ 
son,  i.e.,  Dt,  Dl,  and  Dh>  The  sites  are  listed  in  order  of  longi¬ 
tude  following  a  westwardly  direction.  Values  of  t  were  inter¬ 
polated  to  the  nearest  5  m.  atm-cm  from  the  t  contours  in  Fig. 

I.  Selection  of  the  sites  was  made  initially  from  crude  topo¬ 
graphic  maps.  Exact  site  elevation  values  were  obtained  from 
the  Encyclopedia  Brltannica. 

The  highest  h  value  for  an  inhabited  site  was  found  in 
Cerro  de  Pasco,  Peru,  a  mining  town  of  5,720  people  (1961) 
living  at  an  altitude  of  14,436  ft  or  4.40  km.  The  calculated 
relative  DUV  value  of  1.19  for  Cerro  de  Pasco  was  the  highest 
of  the  19  selected  sites,  even  though  it  was  located  8°  south 
of  the  equator.  The  second  highest  D  value  of  1.17  was  found 
in  Quito,  Ecuador  which  has  the  relatively  modest  altitude  of 
2.85  km  but  lies  on  the  equator.  These  high-altitude  sites  were 
found  to  have  a  relative  DUV  value  approximately  80  percent 
higher  than  Key  West,  Florida. 
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The  simple  model  of  relative  DUV  for  tropical  sites  derived 
in  this  paper  can  be  used  to  tentatively  identify  interesting 
sites  for  the  installation  of  DUV  measuring  instruments.  Site 
selection  would  also  depend  on  ecological  considerations  such 
as  vegetation,  human  habitation,  animal  life,  etc.  Ideally, 
the  selection  of  such  sites  should  not  be  made  until  cloud  cover 
statistics  are  available  for  the  site  candidates.  No  attempt  to 
rank  order  the  19  sites  of  Table  3  by  DUV  was  made  because  of  the 
recognition  that  cloud-cover  statistics,  had  they  been  available 
and  included,  would  almost  certainly  have  altered  such  a  ranking. 

TABLE  3.  RELATIVE  ANNUAL  DUV  DOSE  AT  SELECTED  TROPICAL  SITES 


SITE 

LONGITUDE 

LATITUDE 

t (m. atm-cm ) 

h  ( km ) 

D 

I 

T 

°L 

0 

Entebbe,  Uganda 

32*s°W 

0° 

240 

1  . 

15 

1  . 

00 

1  . 

00 

1 

.07 

1 

.07 

Na i rob 

i ,  Kenya 

37°W 

2°S 

240 

1 . 

68 

1  . 

00 

1 . 

00 

1 

.10 

1 

.10 

Reci fe 

,  Brazil 

35°E 

8°S 

240 

0 

1  . 

00 

0. 

98 

1 

.00 

0 

.98 

Caraca 

s,  Venezuela 

67°E 

1 0°N 

245 

0. 

91 

0. 

98 

0. 

96 

1 

.05 

0 

.99 

La  Paz 

,  Bolivia 

68°E 

1  6°S 

255 

3. 

66 

0. 

93 

0. 

91 

1 

.22 

1 

.03 

Oru  ro , 

Bolivia 

67°E 

1  8°S 

255 

3. 

71 

0. 

93 

0. 

89 

1 

.22 

1 

.01 

Cuzco , 

Peru 

72°E 

1  3*5°  S 

250 

3. 

35 

0. 

95 

0. 

93 

1 

.20 

1 

.06 

Bogota 

,  Columbia 

74°E 

6%°S 

245 

2. 

61 

0. 

98 

0. 

98 

1 

.16 

1 

.11 

Huancayo,  Peru 

75°E 

1 2°S 

250 

3. 

26 

0. 

95 

0. 

95 

1 

.20 

1 

.08 

Cerro 

de  Pasco,  Peru 

76°E 

10°S 

245 

4. 

40 

0. 

98 

0. 

96 

1 

.26 

1 

.19 

Lima, 

Peru 

77°E 

1 2°S 

250 

0. 

15 

0. 

95 

0. 

95 

1 

.01 

0 

.91 

Quito, 

Ecuador 

78>s0E 

0° 

240 

2. 

85 

1 . 

00 

1 . 

00 

1 

.17 

1 

.17 

Panama 

Canal  Zone 

79*s°E 

9°N 

240 

0 

1 . 

00 

0. 

97 

1 

.00 

0 

.97 

Key  We 

st,  Florida 

82°E 

24Ji°N 

275 

0 

0. 

83 

0. 

80 

1 

.00 

0 

.66 

Mauna 

Loa,  Hawaii 

1  56°W 

1  9°N 

260 

4. 

00 

0. 

90 

0. 

87 

1 

.24 

0 

.97 

Townsv 

ille,  Australia 

147°W 

1 9°S 

260 

0 

0. 

90 

0. 

87 

1 

.00 

0 

.78 

Saigon 

,  Vietnam 

1 07°W 

9°S 

255 

0 

0. 

93 

0. 

97 

1 

.00 

0 

.90 

Singapore 

1  04°W 

2°N 

250 

0 

0. 

95 

1 . 

00 

1 

.00 

0 

.95 

Bangal 

ore,  India 

7  7  >a°  W 

1 3°N 

250 

0. 

95 

0. 

95 

0. 

94 

1 

.06 

0 

.95 

III.  MID-LATITUDE  SITES 


A.  DUV  VS  OZONE  THICKNESS  AT  THE  EQUATOR 

To  derive  a  formula  for  relative  DUV  for  mid-latitude  sites, 
it  is  first  necessary  to  modify  Eq.  (4)  for  the  ozone  thickness 
factor  Dt-  The  latter  covered  only  the  tropical  range  of  ozone 
values  240-280  m.  atm-cm.  Since  t  values  can  reach  as  high  as 
400  m.  atm-cm  in  the  Northern  Hemisphere,  it  is  necessary  to  de¬ 
rive  another  formula  for  for  application  to  mid-latitude  sites. 

In  Fig.  17  is  shown  the  annual  relative  DUV  at  the  equator 
as  a  function  of  ozone  thickness  t.  The  solid  line  is  drawn 
through  the  circled  points  which  are  based  on  calculations  of 
A.E.S.  Green  and  T.  Mo  (Ref.  11).  The  relative  DUV  units  in 
Fig.  17  are  based  on  daily  erythema  doses  on  the  15th  day  of 
each  month  In  (Joules/m2 )e  and  differ  from  those  in  Fig.  4 
which  were  based  on  annual  monthly  sums  in  MED. 

A  very  good  parabolic  approximation,  obtained  by  fitting 
the  calculated  DUV  values  at  300,  350,  and  400  m.  atm-cm,  is 
given  by 

Dt  *  3-80  x  10~5  t2  -  3.95  x  10-2  x  +  11.19  .  (22) 

In  order  that  Eq.  (22)  be  compatible  with  Eq.  (4),  the 
equations  are  matched  at  x  =  250.  According  to  Eq.  (4),  Dt  = 
0.9516  at  x  =  250,  whereas,  according  to  Eq.  (22),  Dt  =  3.69. 

Thus,  multiplying  Eq.  (22)  by  0. 9516/3.69  yields  the  normalized 
parabolic  approximation 

Dt  =  9-80  x  10"6  x2  -  1.0186  x  10"2  x  +  2.886  .  (23) 


f 
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FIGURE  17.  Relative  DUV  at  equator  vs  ozone  thickness 


B.  MID-LATITUDE  CORRECTION  FOR  TROPICAL  LATITUDE  FACTOR 

The  tropical  latitude  exponential  factor  DL  in  Eq .  (5) 
must  be  modified  for  mid-latitude  application.  If  DLm  (t,  L) 
denotes  the  mid-latitude  relative  annual  DUV  for  ozone  thick¬ 
ness  t  and  latitude  L,  and  a  (x,  L)  denotes  the  correction 
factor  for  the  value  of  DL  as  given  by  Eq.  (5),  then 

DLm  (T>  L>  ’  [l  -  “  <*’  L>]  dl  <2' 

=  [l  -  a  (T,  L)]  a'3'711  *  10'1‘l2  •  (2; 


In  Fig.  18  are  shown  the  results  of  calculations  of  a  (x,  L) 
based  on  tables  in  Refs.  2  and  11  for  x  values  of  256,  300,  320, 
350,  and  400  m.  atm-cm  and  L  values  of  30°,  35°,  40°,  45°,  50°, 
and  55° •  According  to  Fig.  1,  the  o  values  of  interest  will  lie 


M 


( T,L)  CORRECTION  FACTOR  FOR  Dl 


between  the  two  dashed  lines  in  Fig.  18.  It  is  seen  that  for  a 
given  latitude,  a  is  empirically  found  to  be  a  linear  function 
of  t.  The  slopes  of  the  lines  in  Fig.  18  are  seen  to  increase 
with  latitude.  If  f  (L)  denotes  the  correction  factor  for 
t  =  256  m.  atm-cm,  and  S  (L)  denotes  the  slope  as  a  function  of 
latitude,  then  the  correction  factor  a  (t,  L)  is  given  by 

a  ( t ,  L)  =  f  (L)  +  (t  -  256)  S  (L)  .  (26) 


FIGURE  18.  Correction  factor  for 
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By  fitting  the  values  of  a  (t,  L)  at  t  ■  256  m.  atm-cm  for  lati¬ 
tudes  30°,  45°,  and  55°  to  a  parabola,  it  is  found  that  a  good 
approximation  for  f  (L),  as  illustrated  in  Fig.  19*  is  given  by 

f  (L)  =  9.08  x  10“5  L2  -  5.28  x  10-3  L2  +  7.67  X  10"2  .  (27) 

The  slope  S  (L)  is  found  to  be  linear  in  the  interval 

30°  <  L  <  45°,  as  indicated  in  Fig.  20.  The  function  S  (L)  may 

be  approximated  by  the  equation 

S  (L)  =  g  (L)  +  6  p  (L)  ,  (28) 

where  6=0  for  the  interval  30°  <  L  <  40°  and  6=1  for  L  >  45° 
p  (L)  is  a  slope  correction  factor  for  latitudes  greater  than 
45°.  By  assuming  the  slope  correction  factor  p  (L)  increases 
proportionally  to  the  square  of  (L  -  45°),  the  function  S  (L) 
is  well  approximated  by  the  equation 

S  (L)  =  1.46  x  10~4  +  1.34  x  10-5  (L  -  30°) 

+  6.10  x  10"7  6  (L  -  45 ) 2  .  (29) 

Substituting  Eqs.  (27)  and  (29)  in  Eq.  (26),  the  mid-latitude 
correction  factor  a  (t,  L)  is  therefore  given  by 

a  (t,  L)  -  9.08  x  10-5  L2  -  5-28  x  10-3  L2  +  7.67  x  10-2 

+  (t  -256)  [l.46  x  10-4  +  1.34  x  10-5  (L  -  30°) 

+  6.10  x  10-7  6  (L  -  45°)2]  ,  (30) 

where  6=0  for  the  interval  30°  <  L  <  45°  and  6=1  for 
L  >  45°. 
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C.  CORRECTION  FACTOR  FOR  SEASONAL  OZONE  VARIATION 

If  the  average  ozone  thickness  were  a  constant  value  in¬ 
dependent  of  season,  then  the  substitutions  of  Eq.  (23)  for  the 
ozone  thickness  factor  DT  and  Eqs.  (25)  and  (30)  for  the  lati¬ 
tude  factor  Dl  in  the  tropical  relative  DUV  formula  of  Eq.  (5) 
would  accurately  estimate  the  relative  DUV  for  a  mid-latitude 
site.  However,  the  ozone  thickness  varies  significantly  from 
month  to  month  at  mid-latitude  sites.  This  seasonal  ozone  var¬ 
iation  introduces  a  significant  fractional  error  p  such  that 

Ds  =  D  [l  +  p  (t,  L)]  ,  (31) 

where  D  is  the  relative  annual  DUV  as  adjusted  for  seasonal 
ozone  variation.  It  will  be  shown  below  that  p  is  positive  in 
the  Northern  Hemisphere  and  negative  in  the  Southern  Hemisphere. 
It  will  also  be  shown  that  p  is  a  function  of  the  average  an¬ 
nual  ozone  thickness,  x,  as  well  as  latitude,  L,  in  the  Northern 
Hemisphere,  but  appears  to  be  essentially  independent  of  t  in 
the  Southern  Hemisphere. 

In  the  global  ozone  distribution  atlas  of  Ref.  1,  ozone 
thickness  Is  tabulated  for  each  month  of  the  year  between  July 
1957  and  June  1967,  for  5°  latitude  intervals  in  the  Northern 
and  Southern  Hemispheres,  and  for  20°  longitude  intervals  in 
both  hemispheres.  For  a  given  geographic  location,  let  the  tab¬ 
ulated  entry  x..  denote  the  ozone  thickness  during  the  jth  month 
th 

of  the  i  year.  The  ten-year  average  ozone  thickness  for  the 
jth  month  is  then  given  by 


i=10 


1=1 


and  the  average  annual  ozone  thickness  by 


35 


T 


1 

12 


J-12 


£ 


(33) 


The  average  annual  relative  DUV  dose  D  can  be  accurately 
calculated  by  summing  the  monthly  relative  DUV  dose  using  the 
appropriate  monthly  value  of  ozone  thickness,  ,  i.e., 


D 


s 


3  =  12 

J2  DJ  (V 

3=1 


(34) 


An  approximate  calculation  for  relative  annual  DUV  dose  can  also 
be  obtained  by  using  the  average  annual  ozone  thickness,  t,  from 
Eq.  (33)  for  each  month  of  the  year,  i.e., 

J  =  12 

D  -  ^  Dj  (t)  .  (35) 

J-l 

The  approximate  form  Eq.  (35)  is  the  one  that  was  used  in  the 
tropical  relative  DUV  formula  and  was  based  on  the  tabulated 
monthly  values  given  in  the  tables  of  Refs.  2  and  11.  The  frac¬ 
tional  error  p  that  results  if  seasonal  ozone  variation  is  ne¬ 
glected  and  Eq.  (35)  is  used  for  mid-latitude  sites  is  given 
by 


To  investigate  the  behavior  of  p  with  latitude  and  ozone 
thickness,  calculations  for  D0  and  D  were  made  for  six  sites  in 
the  Northern  Hemisphere  and  five  sites  in  the  Southern  Hemi¬ 
sphere.  The  values  of  longitude  and  latitude,  and  the  calcu¬ 
lated  values  of  t  and  p  for  these  eleven  sites  are  given  in 
Table  4.  In  Fig.  21  the  values  of  p  are  plotted  as  a  function 
of  latitude.  Note  that  in  the  Northern  Hemisphere  the  value  of 
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p  is  positive  and  is  much  greater  at  longitude  100°E  than 
longitude  100°W,  and  increases  with  latitude  at  a  faster  rate. 

In  the  Southern  Hemisphere,  the  value  of  p  is  negative  and  the 
values  for  longitude  60°E  are  very  nearly  the  same  as  for  longi¬ 
tude  60°W. 

TABLE  4.  t  AND  p  VALUES  FOR  ELEVEN  SELECTED  SITES 


Latitude 

Longitude 

T 

P 

25°S 

60°W 

271.1 

-0.015 

45°S 

60°W 

303.7 

-0.033 

55°S 

60°W 

313.7 

-0.053 

45°S 

60°  E 

334.5 

-0.036 

55°S 

CT1 

o 

o 

m 

339.7 

-0.076 

25°N 

100°W 

274.3 

-0.015 

45°N 

100°W 

351.6 

0.030 

55°N 

100°W 

379.4 

0.060 

25°N 

1 00°E 

259.3 

-0.010 

45°N 

1  00°E 

308.9 

0.082 

55°N 

I  00°E 

333.1 

0.164 

Insights  into  the  complexity  of  the  behavior  of  the  seasonal 
ozone  variation  correction  factor  p  can  be  gained  from  compari¬ 
sons  of  the  ten-year  average  monthly  ozone  thickness  Tj  for  the 
eleven  sites  in  Figs.  22a  through  22f.  Observe  the  following: 

1.  In  the  Southern  Hemisphere,  a  minimum  value  of 

occurs  at  mid-latitudes  in  April  and  a  maximum  value 
in  November;  in  the  Northern  Hemisphere,  is  maximum 
in  February  to  March  and  minimum  in  August  to  October. 


FIGURE  22.  Seasonal  variation  of  ten-year  average  ozone 
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2.  The  amplitude  of  the  curve,  i.e.,  the  maximum  value 
minus  the  minimum  value,  increases  with  increasing 
latitude,  is  much  greater  for  the  Northern  than  the 
Southern  Hemisphere  for  a  given  latitude,  and  varies 
more  with  longitude  in  the  Northern  Hemisphere  than 
in  the  Southern. 

3.  The  curves  are  neither  smooth  nor  symmetrical,  indi¬ 
cating  phase  shifts,  e.g.,  in  Pig.  22f,  differences  in 
the  period  for  which  ij  is  greater  than  x,  and  dif¬ 
ferences  in  the  values  of  |xj  -  x|  for  the  maximum  and 
minimum  values  of  Xj  (e.g..  Pig.  22b). 

Despite  the  irregularities  in  the  seasonal  ozone  thickness 
curves,  there  is  also  much  commonality  in  their  characteristics, 
leading  to  fairly  well-defined  seasonal  effects  on  the  relative 
DUV  calculations.  In  Fig.  23,  the  relative  DUV  on  the  15th  day 
of  each  month  is  plotted  for  each  month  of  the  year  with  seasonal 
ozone  variation,  i.e.,  D  (x.^),  and  for  a  constant  annual  mean 
amount  of  ozone,  i.e.,  D  (x),  and  for  each  of  the  five  selected 
sites  in  the  Southern  Hemisphere;  in  Fig.  24  D  ( t ^ ^ )  and  D  (x) 
are  shown  for  each  of  the  six  selected  sites  in  the  Northern 
Hemisphere . 

The  following  conclusions  can  be  drawn  from  the  behavior 
of  the  D  ( Ti j )  and  D  ^  curves  in  Figs.  23  and  24. 

1.  In  the  Southern  Hemisphere,  all  of  the  D  (x^j)  and 

D  (x)  curves  reach  a  minimum  value  in  June.  However, 
and  more  Importantly,  the  D  (x^ )  curves  reach  a  max¬ 
imum  value  in  January,  whereas  the  D  (x)  curves  reach 
a  maximum  value  in  December. 

2.  In  the  Northern  Hemisphere,  the  D  (x^)  are  found  to 
reach  a  maximum  value  in  July,  whereas  D  (x),  of 
course,  always  will  reach  a  maximum  value  in  June. 

3.  The  symmetry  of  the  D  (x)  curves  is  not  to  be  found 

in  the  D  (x^)  curves,  particularly  in  the  Northern 
Hemisphere.  j.n 
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FIGURE  23.  Comparison  of  dally  DUV  using  seasonal  ozone  variation 
and  using  constant  mean  value  of  ozone  for  five  sftes 
in  the  Southern  Hemisphere 
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In  the  Southern  Hemisphere,  from  January  to  June, 
D  (fjj)  >  D  (x);  from  July  to  December,  D  (x^)  < 
D  (t). 


5.  In  the  Northern  Hemisphere,  from  January  to  May  or 
June,  D  (x^)  K  D  ( t ) ;  from  June  or  July  to  December, 

D  (x^)  >  (x). 

6.  The  values  of  p  are  not  very  large,  because  the  dif¬ 
ference  in  the  D  (t^j )  and  D  (x)  curves  during  the 
first  half  of  the  year  is,  to  a  large  extent,  negated 
by  the  difference  in  the  two  curves  of  opposite  sign 
during  the  second  half  of  the  year. 

For  the  Southern  Hemisphere,  it  is  possible  to  approximate 
p  as  a  function  of  latitude  with  a  single  equation  (Fig.  21). 

By  using  the  mean  value  of  p  for  longitudes  60°E  and  60°W  at 
latitudes  45°S  and  55°S,  a  parabolic  fit  of  p  at  latitudes 
25°S,  45°S,  and  55°S  results  in  a  seasonal  ozone  correction 
factor  ps  for  the  Southern  Hemisphere  of 


po  (L)  =  -  6.75  x  10~5  L2  +  3.75  x  10  3  L  -  6.66  x  10-2, 
s 


(37) 


As  can  be  seen  from  inspecting  Fig.  1,  longitude  60°W  lies 
in  a  region  of  low  values  of  x  at  mid- latitudes ,  whereas  longi¬ 
tude  60°E  lies  in  a  region  of  relatively  high  values.  It  would 

therefore  appear  to  be  the  case  that  p_  is  essentially  independ- 

s 

ent  of  x  or  the  seasonal  ozone  fluctuations  that  tend  to  be 
highly  correlated  with  x. 

In  the  Northern  Hemisphere,  the  seasonal  ozone  correction 
factor  pn  is  strongly  dependent  on  x  as  well  as  L  and  so  a  more 
complex  formulation  Is  required. 

The  behavior  of  x  in  the  Northern  Hemisphere  with  a  change 
in  longitude  is  shown  in  Fig.  25.  Note  that  longitude  100°E, 
which  at  mid-latitude  represents  a  region  of  minimum  amount  of 
ozone  and  seasonal  ozone  fluctuations  is,  surprisingly,  a  region 
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FIGURE  25.  Mean  annual  ozone  thickness  vs  latitude  for 
longitudes  100°E  and  100°W  In  the  Northern 
Hemisphere 
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of  maximum  values  of  the  seasonal  ozone  correction  factor  p 
(Fig.  21),  whereas  longitude  100°W,  which  represents  a  region 
of  maximum  amount  of  ozone  and  relatively  large  seasonal  ozone 
fluctuations,  is  a  region  of  minimum  values  of  p.  Let  x  (L) 
represent  the  ozone  thickness  as  a  function  of  latitude  for 
longitude  100°E.  From  Fig.  25  it  is  seen  that  xq  (L)  can  be 
well  approximated  by  the  linear  equation 


xQ  (L)  =  197.8  +  2.46  L  .  (38) 

The  equation  pQ  (L)  for  p  at  longitude  100°E  in  Fig.  21  can  be 
approximated  by  fitting  the  values  of  p  at  latitudes  25°,  45°, 
and  55°.  The  result  is 

Po  (L)  =  0.010  +  1.20  x  10-11  L2  -  3.80  x  10~3  L  .  (39) 

The  equation  for  p  at  a  longitude  between  100°E  and  100°W, 
and  hence  a  value  of  x  between  the  two  curves  in  Fig.  25,  can 
be  approximated  by  the  equation 


p  (x,  L)  =  p0  (L)  -  [x  -  xQ  (L)]  q  (L)  (40) 


where  q  (L)  is  the  difference  in  the  p  curves  for  longitudes 
100°E  and  100°W  (Fig.  21)  divided  by  the  difference  in  the 
corresponding  x  values  (Fig.  25).  Thus, 


p  (L)  -  p..  (L) 

-  ■  T°-a)  -  r0  (l)  • 


(Hi) 


where  p.^  (L)  and  x^  (L)  are  the  seasonal  ozone  correction  fac¬ 
tor  and  ozone  thickness,  respectively,  for  longitude  100°W. 

By  fitting  a  parabola  to  the  latitude  points  25°,  45°,  and  55°, 
it  can  be  shown  that 


q  (L)  -  1.424  x  10"3  +  1.955  x  10“6  L2  -  9.25  x  10~5  L.  (42) 


45 


Substituting  Eqs.  38,  39,  and  M2  in  Eq.  MO  results  in  the 
following  approximate  formula  for  p  in  the  Northern  Hemisphere 
as  a  function  of  L  and  t. 

Pn  (t,  L)  =  0.010  +  1.20  x  10"21  L2  -  3.80  x  10'3  L  (M3) 

-  (x  -  197.8  -  2.M6  L)  (1.M2M  x  10-3  +  1.955  x  10~6  L2 

-  9.25  x  10"5  L) 

D.  FORMULA  FOR  RELATIVE  ANNUAL  DUV  AT  MID-LATITUDE  SITES 

Substituting  Eq.  (23)  for  the  first  term  in  Eq.  (21),  mul¬ 
tiplying  by  the  latitude  correction  [1  -  a  (x,  L)]  and  the  sea¬ 
sonal  ozone  correction  factor  p  gives,  for  the  relative  annual 
DUV  at  mid-latitude  sites,  the  formula 

D  *  (9-80  x  10-6  x2  -  1.0186  x  10"2  x  +  2.886) 

-M2 

x  [1  -  a  (x,  L)]  e“3’ x  10  L  [1  +  p  (x,  L)]  (1  +  0.06  h) 
x  (1  -  0.50C)  (1  +  0.50  A)  [1  -  0.093  (S  -  1)]  (MM) 

where  a  (x,  L)  is  given  in  Eq.  (30),  and  p  (x,  L)  by  Eq.  (37) 
for  the  Southern  Hemisphere  and  by  Eq.  (M3)  for  the  Northern 
Hemisphere.  Other  symbols  are  as  defined  on  p.  26. 
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IV.  CONCLUSIONS 


It  is  shown  in  this  paper  that  it  is  possible  to  derive 
and  apply,  to  the  extent  that  available  ozone  and  meteorologi¬ 
cal  and  geographical  data  are  available,  a  formula  for  a  rela¬ 
tively  quick  determination  of  annual  relative  DUV  dose  for  trop¬ 
ical  and  mid-latitude  sites.  The  input  parameters  required  are 
average  annual  amount  of  ozone,  latitude,  altitude,  average 
cloud  amount,  ground  albedo,  and  amount  of  aerosols.  Ozone, 
latitude,  and  altitude  information  is  readily  available. 

Cloudiness  information  to  the  accuracy  required  is,  unfortunately, 
not  readily  available  on  a  worldwide  basis.  Input  data  on  ground 
albedo  and  aerosol  content  are  also  not  readily  available,  but, 
in  general,  can  be  expected  to  play  a  less  significant  role  than 
the  other  four  parameters. 

The  derived  formulas  for  tropical  and  mid-latitude  sites 
can  be  used  to  provide  a  fundamental  input  parameter,  relative 
annual  DUV  dose,  in  models  designed  to  investigate  the  effects 
of  solar  ultraviolet  radiation  and  its  possible  increase  result¬ 
ing  from  stratospheric  ozone  depletion,  on  ecological  systems 
on  land  and  in  the  oceans  of  the  world,  and  on  the  incidence  of 
skin  cancer  in  white  Caucasian  populations. 

Because  of  the  sparsity  of  the  data  available  and  possible 
interdependency  effects  in  the  six  multiplicative  factors  de¬ 
rived,  the  formula  remains  to  be  validated  in  a  general  sense. 
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